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Abstract. We apply the slave-boson approach of Kotliar and Ruckenstein to the two-band Hubbard model 
with an Ising like Hund's rule coupling and bands of different widths. On the mean-field level of this 
approach we investigate the Mott transition and observe both separate and joint transitions of the two 
bands depending on the choice of the inter- and intraorbital Coulomb interaction parameters. The mean- 
field calculations allow for a simple physical interpretation and can confirm several aspects of previous 
work. Beside the case of two individually half- filled bands we also examine what happens if the original 
metallic bands possess fractional filling either due to finite doping or due to a crystal field which relatively 
shifts the atomic energy levels of the two orbitals. For appropriate values of the interaction and of the 
crystal field we can observe a a band insulating state and a ferromagnetic metal. 



PACS. 7f .30.-|-h Metal-insulator transitions and other electronic transitions - 
and alloys - 11.27.+& Strongly correlated electron systems; heavy fermions 
models (Hubbard model, etc.) 



7f .20.Be Transition metals 
- 71.fO.Fd Lattice fermion 



1 Introduction 

The Mott transition in multiorbital systems with several 
bands gives rise to complex and intriguing physics. Multi- 
band systems occur naturally in rare earth intermetallic 
compounds and in systems involving transition metals. In 
the former extended conduction electrons and almost lo- 
calized /-electrons couple through local hybridization and 
give rise to Kondo and heavy Fermion physics. In transi- 
tion metal oxides, chalcogenides etc. several partially filled 
d-orbitals are the origin of rather similar electron bands of 
different but comparable width. Here the question arises 
how the interaction among these orbitals. Coulomb repul- 
sion and Hund's rule coupling, influences the transition 
to partially or fully localized degrees of freedom. What is 
the nature of the Mott transition that occurs as the mag- 
nitude of the interactions is increased gradually? In this 
paper we will be concerned with these questions. 

The two-band Hubbard model with bands of different 
widths is the simplest model that captures all the relevant 
aspects of the Mott transition in multiorbital systems. In 
recent years this model was investigated by several authors 
J., 6, 7, 8, 9, lOj li] mainly in the framework of 
dynamical mean-field theory (DMFT) and using different 
methods to solve the local impurity problem. These calcu- 
lations have led to the following understanding of the Mott 
transition at half filling. Depending on the exact choice of 
the intra- and interorbital interaction parameters, one can 
observe a sequence of individual Mott transitions in each 
band or a joint transition involving both bands simulta- 



neously. The existence of a separate transition, usually re- 
ferred to as "orbital-selective Mott transition" (OSMT), 
implies an intermediate phase between the metal and the 
Mott insulator where only the narrow band is insulating 
whereas the wide band still has metallic properties. Fur- 
thermore, the stability of this intermediate phase strongly 
depends on how the Hund's rule coupling is taken into 
account. 

Early studies of the Mott transition in multiorbital sys- 
tems made by Anisimov et al. Liebsch 0,0,0,01 and 
by Koga and coworkers 0, 0, lorand more recent DMFT 
calculations of de' Medici et al. Ferrero et al. 0, Arita 
et al. , and Knecht et al. showed that different im- 
purity solvers capture different aspects of the Mott tran- 
sition and can partially lead to different conclusions. It 
is therefore desirable to investigate the properties of this 
Mott transition within a more analytical theory. We a pply 
the slave-boson approach of Kotliar and Ruckenstein [13 
on the mean-field level, discuss and confirm several aspects 
of previous work. Our calculations give reasonable results 
in a wide range of parameters and allow in a natural way 
for a simple physical interpretation. Beside the case of two 
individually half-filled bands we examine what happens if 
the original metallic bands possess fractional filling either 
due to finite doping or due to a crystal field which shifts 
the atomic energy levels of the two orbitals relative to 
each other. In both cases we can observe an OSMT. Due 
to the crystal field splitting also a ferromagnetic and a 
band insulating phase appear in the phase diagram. 
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The paper is organized as follows. In Sec. [51 we intro- 
duce the model. In the rather technical Sec.|21we apply the 
slave-boson formalism and its mean-field approximation. 
The results concerning the Mott transition are presented 
and discussed in Sec. 0] for various choices of the interac- 
tion parameters at half filling, in the presence of a crystal 
field and for finite doping. Conclusions and a comparison 
with previous results are found in Sec. El 



2 Model 

We consider the following two-band Hubbard Hamiltonian 

^ = EE*?4.c.- + F (1) 

aa (ij) 

with 

V = U ^nialflial +'^iU' - J5„„>)hii„ni2a' ■ (2) 

As usual cl^.^. [ciaa) crcatcs (annihilates) an electron with 
spin a =T,i and band index a = 1,2 at the site i and 
niaa — cl^a-'^iaa is the Corresponding occupation num- 
ber operator. The hopping integral for the orbital a is 

denoted by t|"^ . We assume vanishing interorbital hy- 
bridization and that the hopping integrals have different 
values for the different orbitals, i.e. that the tight-binding 
bands have different bandwidths. The intraband (inter- 
band) Coulomb repulsion is denoted by U ([/') and the 
Hund's rule coupling by J. In two-band Hubbard models 
additional spin-flip and pair-hopping terms are usually in- 
cluded in the Hund's rule coupling. As shortly discussed in 
the next section, these terms pose problems in the slave- 
boson formalism and we therefore concentrate on the Ising 
like Hund's rule coupling in Note however that these 
terms are not included in Quantum Monte Carlo (QMC) 
calculations either [iJ, |^ . For a spherically symmetric 
screened Coulomb interaction the positive interaction pa- 
rameters are related by U' — U ~ 2J 14]. The relevant 
parameter regime is therefore U > U' where the intraor- 
bital repulsion is bigger than the interorbital. 



3 Slave-boson formulation of the two-band 
Hubbard model 

3.1 Slave-boson model 

The treatment of on-site interactions with slave bosons is a 
well established method in different fields of strongly cor- 
related electron systems. Kotliar and Ruckenstein [l^ in- 
troduced this approach for the (one-band) Hubbard model. 
By a new functional-integral representation of the parti- 
tion function they were able to effectively map the fermionic 
action on a bosonic action with local constraints. The 
simplest saddle-point approximation of their approach re- 
produces the results of the Gutzwiller approximation |l(l 



Table 1. The atomic states in the original model, their corre- 
sponding slave-boson states as well as the labeling of the mean 
fields. The site index is suppressed, a = 1,2, and a =J. (t) if 

'^=T (i). 
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llTl Il8| | . The slave-boson approach leads to a novel mean- 
field theory which is especially useful for examining the 
Mott transition. As the Gutzwiller approximation, the 
slave-boson mean-field theory is closely related to Lan- 
dau's Fermi liquid theory 0| since the slave bosons keep 
track of the other electrons by measuring the electron oc- 
cupancy at each atom which leads to a renormalization of 
the hopping amplitude and thus to a change of the effec- 
tive mass. Let us first look at one particular lattice site. 
The atomic Hilbert space is 16-dimensional and spanned 
by the local occupation number basis listed in the first 
column of Tab. and sketched in Fig. The essence of 
the slave-boson approach of Kotliar and Ruckenstein is to 
map the original fermionic model to a mixed fermionic- 
bosonic model with local constraints by introducing for 
each atomic configuration an auxiliary boson 

where a = 1,2, a —1, |. The labeling of the boson opera- 
tors is sketched in Fig. ^ In the following we denote the 
fermionic annihilation (creation) operators in the slave- 
boson model by f^^^ to distinguish them from c^^^^. defined 
in the purely fermionic model. In the extended model, the 
creation of a general slave-boson state is realized by act- 
ing with the bosonic creation operators © and the new 
fermionic operators on the vacuum |vac). The states which 
correspond to the physical atomic states of the original 
model are listed in the second column of Tab. ^ 

The introduction of the bosonic degrees of freedom 
leads to unphysical states which are eliminated by local 
constraints. Summing up all boson occupancy operators, 
we define 

acr 

+ E + E dl.'d^aa- + ///,. (4) 
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Fig. 1. The sixteen atomic configurations of the two-band Hubbard model and the corresponding slave bosons. 



Furthermore we define the operators 

a' 

+ hlM<^ + Y.^A2a'h^2a' + flh. (5) 
a' 

a' 

+ hLh^2a + J2 ^L'^'l-' + (6) 



where 



The physical subspace is given by the local constraints 



(7) 
(8) 



These constraints ensure that the slave-boson states listed 
in the second column of Tab. ^ form a complete set in the 
physical local Hilbertspace of the slave-boson model. The 
first relation ((JJ) represents the completeness of the bo- 
son operators, i.e., the sixteen states with one boson form 
a complete set in the local physical Hilbertspace of the 
bosons. The operators Qiao- count the number of bosons 
that correspond to local configurations having an electron 
with spin a in the orbital a. Therefore, we have to ensure 
with the constraint © that in the physical subspace the 
operators Qiaa are identical to the operators fj^^fiaa- 

Using these constraints, the interaction term becomes 
quadratic in the boson operators 

= E E + + - E ^\o.ah^.a 

+{u' dLd,,^ + u'J2 dl-Aas 

+2{U + 2U' -J)flhy (9) 

The attempt to include the spin-flip and pair-hopping 
term in a similar way fails due to quartic fermion terms 
in V^^ or in the constraints. In this case additional ap- 
proximations are required Whereas the interaction 
term has become much simpler the new formulation of 
our model implies that the destruction or creation of a 
physical fermion has to be accompanied by slave bosons. 



ft 

J taa Ida ' 



(10) 



Zi2a = e\pi2a + p\2sSi2 + p\i„di„„ + Pagdiacr 

+ 4lht2a + 4„ghaa + 4a9hrla + 429 f^' 



The "z-operators" keep track of the environment (bosons) 
during hopping processes The choice of the "z-operators" 
is not unique. In fact, as long as the constraints are ful- 
filled exactly, the "z-operators" can be modified by any 
operator which is the identity operator when restricted to 
the physical subspace ^15, 20] . However, the mean-field re- 
sults depend on the choice of these operators. The choice 
of Kotliar and Ruckenstein, that we take here HlOfl. repro- 
duces correctly the noninteracting case in the mean-field 
approximation 15]. The slave-boson Hamiltonian is then 
given by 



ijaa 



(11) 



and is fully equivalent to the original Hamiltonian pro- 
vided the local constraints (|7||H1) are handled exactly. They 
can be imposed by site dependent Lagrange multipliers 
and X?„„. 



3.2 Mean-field approximation 

The simplest saddle-point approximation of the grand canon- 
ical partition function Z = Tr[e^'^'^ "''^ ^V] is equiva- 
lent to a mean-field approximation where the Bose fields 
and Lagrange multipliers are treated as static and homo- 
geneous fields. Thus, this approximation consists essen- 
tially in replacing the creation and annihilation operators 
of the slave bosons by site independent c-numbcrs which 
can be chosen to be real. The mean fields are listed in the 
fourth column of Tab. ^ In this approximation, the con- 
straints are fulfilled only on average and the square of the 
mean fields can be interpreted as the probability of find- 
ing the corresponding local configuration at a particular 
site. The mean-field Hamiltonian with included averaged 
constraints can be diagonalized and yields at T = the 
variational ground-state energy (per site) 

acr aa 

(12) 
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By introducing the effective chemical potential fiaa = M ~ 
A^g. in band a for fermions with spin cr, the average kinetic 
energy per site and band reads 



2.5 



(13) 



Similarly, the density in band a is 



depa{e). (14) 



The chemical potential [iaa has to be determined from 
Eq. (|14|l for a given density. If both bands are half-filled 
separately one finds for example \Xaa = 0. The bare den- 
sity of state (DOS) per spin in the band a is denoted by 
Pa and the mean-field Coulomb energy per site is 



\ <j<j' 



4/^ 



2/2 



(15) 



The band-renormalization factor qaa = can be related 
to the effective mass of quasiparticles of Landau's Fermi 
liquid theory. For quasiparticles in the band a with spin 
cr we have q~l = TO*g./m fl3,l2lj. The vanishing of qaa 
therefore indicates the transition to a localized state. The 
mean fields and the Lagrange multipliers are determined 
by the stationary point of Eq which is a saddle point but 
not a minimum. With the help of the averaged constraints 
and of Eq. H14I) we reduce the number of independent vari- 
ables. The stationary point of the variational ground-state 
energy per site 



Eg — ^ Qaa^aa 



MF 



(16) 



then becomes a true minimum and can be found numeri- 
cally in a rather simple way. 



4 Mott transition in the two-band model 
within slave-boson theory 

In this section we present the results concerning the Mott 
transition obtained by numerically minimizing Eq. (|16(l . 
Unless otherwise stated we always assume a paramagnetic 
ground state, i.e. Ua^ = Uai, and a particle-hole sym- 
metric bare DOS. Consequently, the following conditions 
are satisfied: e — f , Pa = Paa — haa', s = si — S2, 
do = d^i = dii and di = o?|| ~ d^. This greatly reduces 
the computational effort. 

If both bands are separately half-filled, the results of 
the mean-field calculations do not depend on the exact 
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Fig. 2. Phase diagram at half filling for U > U' , J = and 
D1/D2 = 1/2. Three diflferent phases can be distinguished: a 
paramagnetic metal (PM), a Mott insulator (MI) and in be- 
tween an orbital selective Mott insulator (OSMI). Two second 
order lines (dashed) merge at Pt to a single first order line 
(solid) which ends in a critical second order point Pc- 



choice of the bare DOS, as long as it is particle-hole sym- 
metric. Away from half filling, the mean-field results slightly 
depend on the exact choice. For simplicity, we choose through- 
out this section for both bands a rectangular DOS and 
denote its half- width by Da- The narrow band is always 
referred to as band 1 and the wider as band 2. Unless 
otherwise stated we choose the ratio of the bandwidths 
such that^ D1/D2 — 1/2. Energy is measured in units of 
the bandwidth of band 2, i.e. 2D2 — 1. We restrict to 
the relevant parameter regime U > U' , where the intraor- 
bital repulsion is bigger than the interorbital. In Sec. 14.11 
we discuss the U-U' phase diagram for J = at half fill- 
ing. In Sec. 14.21 we focus on the dependence of the Mott 
transition on the ratio D1/D2 for vanishing Hund's rule 
coupling and U = U'. In Sec. 14. 31 we impose the condition 
U' = U — 2J. In this case we also examine the effect of a 
crystal field and the influence of finite doping. 



4.1 U > U' andJ = 

The phase diagram at half filling is displayed in Fig. [21 
We can distinguish three different phases: a metallic state 
(PM), a Mott-insulating state (MI) and an intermediate 
state (OSMI) induced by the OSMT where the localized 
band 1 coexists with the metallic band 2. For vanishing 
Hund's rule coupling, our calculations suggest that the 
OSMI phase is bounded by two second order lines (dashed 
lines in Fig.|2). They merge to a single first order line at 
Pt which ends in a critical point Pc- The occurrence of 
this orbital selective Mott insulator is not surprising since 
we have neglected local spin-spin interactions. Therefore, 
if the first band is in a Mott-insulating state, the electrons 
in the second band only feel a uniform charge background 

^ Originally, this choice was motivated by the fact that in 
Ca2_2,Sri,Ru04 the d^ji^^andwidth is approximately twice the 
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Fig. 3. Mean fields of the slave bosons and band- 
renormalization factors (a) for U' — J — and (b) for 
U = U' J = 0. 



(arising from the localized electrons in orbital 1) and the 
Mott transition in the broader band occurs at the critical 
interaction strength Uc2 — 8|£2l — 4Z52 = 2 which is the 
value of one independent band with bandwidth 2D2 = 1 
\lK IT^ . To understand the behavior of the system for 
general values < J7' < J7 it is instructive to consider 
first the following two limiting cases: 

i) U' = . In this case the two bands are independent 
and the critical interaction strength of the Mott transi- 
tion is proportional to the bandwidth Da- In Fig. 3(a) we 
have plotted the mean fields of the slave bosons and the 
band-renormalization factors. In the noninteracting sys- 
tem, U — 0, all configurations are equal likely: e = s — 
Pi = P2 = d — 1/4 where d = do = di. The vanish- 
ing of qi is accompanied by the vanishing of e, s and p2, 
whereas (72 becomes simultaneously zero with pi. In the 
Mott-insulating phase {U > Uc2) we find at each lattice 
site one of the possible four atomic configurations rep- 
resented by the boson d and therefore d reaches 1/2 at 

U = Uc2- 

a) U = U' . For this choice the interaction Hamilto- 
nian has an enlarged symmetry with six degenerate two- 
electron configurations shown in Fig. ^ four spin config- 
urations with one electron in each orbital (represented by 
d) and two configurations with both electrons in one of the 
two orbitals (represented by s). This higher symmetry is 



4f 



d s 

Fig. 4. For U = U' and J = Q there are six degenerate two- 
electron on-site configurations. They are represented by the 
mean fields d and s 



due to the fact that the Coulomb energy of a local con- 
figuration depends only on the total charge on the atom. 
The additional symmetry in orbital and spin degrees of 
freedom enlarges the phase space for charge fluctuations 
and leads to a stabilization of the metallic phase . In 



Fig. [3(b) the mean fields of the slave bosons as well as the 
band-renormalization factors are plotted as a function of 
U . There is a joint Mott transition at the critical inter- 
action strength Uc- Because of orbital fluctuations and in 
contrast to the case U' <U not only conflgurations repre- 
sented by d but also conflgurations represented by s have 
a finite probability at Uc- Despite the high symmetry, it 
surprisingly turns out that d ^ s. The relative strength of 
the mean fields in Fig. 3(b) can be understood as follows. 
The high Coulomb energy of a fully occupied local conflgu- 
ration disfavors most strongly the mean field e. The effect 
of the intraband Coulomb interaction is stronger in the 
narrow band and therefore pi > p2 since pi favors local- 
ized behavior in the band 1 and itinerant behavior in the 
band 2. Furthermore, the z-factors can be approximated 
by zi « 2pis + Ap2d and Z2 ~ 2p2S + '^Pid for high values 
of U . In order to optimize the hopping in the wide band, 
d is slightly increased compared to s. Note that above Uc 
the ratio of d and s is not determined at zero temperature. 

The behavior of the system for general values Q <U' < 
U is mostly determined by the physics of the above dis- 
cussed two special choices of parameters. 



4.2 Di < D2, U = U' andJ = 

Within our mean-field calculation the existence of a joint 
transition for U — U' and J — depends on the ratio of 
the bandwidths Di/D2. It turns out that for ratios below 
a critical value {Di/D2)c the mean- field calculations sug- 
gest an OSMT even for U — U' . Thus, we recover exactly 
the same results as Ferrero et al. ^ using the Gutzwiller 
approximation and de' Medici et al. within their slave- 
spin approximation. The critical ratio where an OSMI oc- 
curs can be calculated analytically within the Gutzwiller 
(or slave-boson) approximation and is ^Qj {Di/D2)c = 1/5. 
At first sight, the existence of such a critical ratio seems to 
be in contradiction with the symmetry argument given by 
Koga et al. It states that for vanishing J and U' = U 
the Mott-Hubbard gap in both bands closes at the same 
critical interaction strength, independent of D1/D2. How- 
ever, it does not exclude the transition into an intermedi- 
ate phase, where the "localized" band is not fully gapped 
J]. Indeed, DMFT calculations for < D2 HQ show 
clearly that the "localized" band is not fully gapped but 
has spectral weight down to arbitrarily low energies. This 
subtle aspect is not captured by the Gutzwiller approxi- 
mation and related mean-fleld theories. 

There is however the possibility that the OSMI phase 
is replaced by an instability not considered so far. Our 
mean-field calculation as well as earlier DMFT calcula- 
tions did not take into account a possible enlargement of 
the unit cell. Below a critical temperature Tpf one usually 
finds antiferromagnetic long-range order in the Mott in- 
sulating phase depending on the topology of the lattice. 
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Interestingly for U' — U and J = spin and orbital de- 
grees are relevant and it is possible that the OSMI phase 
is unstable against an orbitally ordered phase. Whereas 
for U' < U there is no tendency toward such an instabil- 
ity it cannot be excluded a priori for U' = U. We discuss 
now the mechanism which may drive orbital order as an 
another way to double the unit cell. 

Let us look at the extreme limit Di U ~ U' D2 
and assume that the lattice is bipartite. In an adiabatic 
approximation the narrow band is localized and fully dom- 
inated by the intraorbital Coulomb repulsion [Di/U ~ 0) 
whereas the intraorbital Coulomb interaction in the wide 
band has negligible effect {U/D2 « 0). We look at the 
following two limiting cases for the static configuration of 
the localized (narrow) band and their implications to the 
electronic properties of the wide band: 

(a) Homogeneous charge distribution with exactly one elec- 
tron per orbital. 

(b) Staggered charge distribution with doubly occupied or- 
bitals on one sublattice and empty orbitals on the other 
sublattice. 

In the case (a) , the homogeneous charge background con- 
tributes an amount U per site to the total energy. In the 
case (b) , the doubling of the unit cell and the induced re- 
arrangement of electrons in the wide band opens a gap. In 
this way the interorbital interaction is reduced (< U). For 
perfect nesting the wide band is fully gapped and shows 
insulating behavior. On the other hand, the doubly occu- 
pied orbitals in the localized phase cost a fixed amount 
U /2 per site. Thus, there is a competition between these 
two effects which can favor an orbitally ordered phase in 
a certain parameter range. 

In summary, the Mott transition for U' = U is a. subtle 
issue due to the aspect of possible orbital order and we 
suggest that this plays a relevant role for the case Di ^ 
I?2- Analogous to the spin ordering in the Mott insulator, 
the stability of such a phase depends also on details of the 
band structures. Taking into account the possibility of a 
doubling of the unit cell is an interesting topic for further 
investigations to be reported elsewhere. 



4.3 U > U' and U' = U - 23 

We turn back to a given ratio D1/D2 = 1/2. From now on 
we adopt the relation U = U' + 2UJ which is usually used 
in the discussion of the Mott transition in the two-band 
Hubbard model. This relation is valid for a rotationally 
symmetric (screened) Coulomb interaction. 



4.3.1 Mott transition at half filling 

The phase diagram is shown in Fig. [51 Again we can ob- 
serve the OSMI phase, but it is limited to a tiny parameter 
regime. In general the Mott transition is shifted to smaller 
values of U since the Ising like Hund's rule coupling favors 
localized configurations with parallel spins. For the same 
reason the Mott transition in the second band is closer 
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Fig. 5. Phase diagram at half filling for U = U' + 2J and 
D1/D2 ~ 1/2. Three different phases can be distinguished: a 
paramagnetic metal (PM), a Mott insulator (MI) and in be- 
tween an orbital selective Mott insulator (OSMI). At Pt the 
second order line (dashed) meets the first order line (solid) 
which ends in a continuous critical point P^. 



to the one in the first band because, in contrast to the 
case J = discussed in Sec. 14.11 the electrons in the sec- 
ond band not only feel a uniform charge background but 
also a localized spin at each lattice site after the gap for 
charge excitations in the first band has opened. For small 
values of J the physics for U' = U becomes important. 
The dashed line in Fig.|5l which separates PM-OSMI, is a 
second order line whereas the solid line, which separates 
OSMI-MI and PM-MI, is a first order line^. They merge 
at Pt . The first order line ends in a second order transition 
point Pe (Fig.|3(b)i. 



To illustrate the occurrence of the OSMI and the first 
order transition line we show in Fig. the slave-boson 
mean fields and the band-renormalization factors for a 
fixed ratio U'/U =1/2 and J/U = 1/4. We clearly see 
that there is a sequence of individual transitions and that 
(72 jumps at U ~ 0.89 from a finite value to zero. The dis- 
continuity is also observed in the mean fields do, and 
pi. We computed the ground-state energy as a function 
of pi for different values of U, where U' and J have the 
same ratio as above. This is shown in Fig.d Note that pi 
represents the probability to find at a particular site one 
electron in orbital 1 and either no electron or two electrons 
in orbital 2 and serves therefore as the order parameter 
for the Mott transition in the second band. At [/ « 0.89 
the metallic solution pi 0.2 and the Mott-insulator so- 
lution pi = are degenerate. This results in a first order 
transition and a finite jump in pi and consequently also 
in (72 (Fig. EJ. 



^ Including spin-flip and pair-hopping terms in the Hund's 
rule coupling Koga et al. 4] reported two successive second- 
order transitions at T = 0. For T > Liebsch Q identified a 
sequence of two first-order transitions for the same model. If 
spin-flip and pair-hopping terms are omitted he found a first- 
order transition followed by a continuous transition. 
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Fig. 7. Ground-state energy as a function of pi for different 
values of U, where U'/U =1/2 and J/U = 1/4. pi is the 
order parameter for the Mott transition in the second band. 
At [/ ~ 0.89 the metallic solution pi « 0.2 and the insulating 
solution pi = are degenerate. This results in a first order 
transition. 



4.3.2 Effect of a crystal field 

Until now we have assumed that the two bands are both 
centered symmetrically around the Fermi energy. What 
happens if a crystal field splits the atomic energy level for 
the two different orbitals? Let us assume that the overall 
system is still half-filled. We introduce an external field 
?7^^(nii — •hi2) in the Hamilton operator which splits 
the atomic energy levels by 21]. Particle-hole symmetry al- 
lows to concentrate on ry > 0. In the noninteracting case, 
this leads to a relative shift of the tight-binding bands by 
2r] and if this shift is bigger than Di + D2 the lower band 
is totally filled whereas the upper band is empty. In this 
case the system is a band insulator. How does this band 
insulator evolve when we turn on the Coulomb interac- 
tion? Is there a transition from the band insulator to the 
Mott insulator, or can we observe a new phase in between? 

To answer these questions we investigate the effect of 
the crystal field on the mean-field level of the slave-boson 
approach. In contrast to the case 77 = we also keep the 
spin dependence of the mean fields so as to detect a pos- 
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Fig. 8. Phase diagram in the presence of an external field 
77 for U'/U = 1/2 and J/U = 1/4. The OSMI is limited to 
a small parameter regime as shown in the inset. The crystal 
field introduces two new phases: a band insulator (BI) and a 
ferromagnetic metal (FM). 



sible ferromagnetically ordered state. The external field 
leads to an additional term in the variational ground-state 
energy Eq ^ + vi^i - ^2). 

Let us first discuss the phase diagram for a finite Hund's 
rule coupling. To be specific we fix U'/U =1/2 and J/U = 
1/4. The result of the minimization of the ground-state en- 
ergy for different values of U and 77 is shown in Fig. |S1 In 
our slave-boson approach the OSMI phase is restricted to 
a tiny parameter regime and only present for small values 
of rj as shown in the inset of Fig. |S1 The transition in the 
first band (dashed) is a second order line which merges 
the first order transition line (solid) at Pf. If the crystal 
field is strong enough, the system is in a band insulating 
state (BI), i.e. one energy band is totally filled whereas 
the other is empty. For rj > this state is characterized 
by the mean field S2 = 1 • For U ~ the transition PM-BI 
is second-order and happens at r/ = {Di + D2)/2 = 0.375. 
For a finite U the transition is first-order since the charge 
abruptly jumps from n2 < 1 to ?t.2 = 1. For very strong 
values of U and rj there is a competition between the BI 
and the MI phase. The boundary is given by comparing 
the energy of the BI phase, U — 2r], with the energy U' ~ J 
of the MI phase and yields U = 877/8 for the above given 
ratio of the interaction parameters. The most interesting 
region of Fig. |H| lies between these limiting cases where a 
ferromagnetically ordered metal (FM) is observed. This 
state is twofold degenerated and triggered by the finite 
Hund's rule coupling. Within our approximation we al- 
ways find maximal spin polarization which is character- 
ized by a finite value of the mean fields {s2, /iio-, rfo-o-,P2o-}- 
We can get an idea of the physical mechanism by fixing 
77 = 0.5 and increasing U continuously starting at t/ = 0. 
The evolution of the charge is shown in Fig. |5| At the 
beginning the Coulomb repulsion is too weak to put elec- 
trons in the upper band and 7i2| = = 1. At a critical 
interaction strength it is energetically favorable to popu- 
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solution (PM) there is a transition to a band insulating phase 
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late the upper band by a few electrons of the same spin 
species and ni| jumps from to a finite value. In the lower 
band simultaneously jumps from 1 to a value < 1. 
The Pauli principle excludes doubly occupied orbitals in 
the upper band which reduces the Coulomb energy and 
leads to a ferromagnetic order. In addition, J couples the 
spin between upper and lower band and we find the same 
magnetization in both bands: toi = rn2 = to/2. Note that 
the critical interaction for the transition to the FM phase 
depends on the exact choice of the bare DOS. Increasing 
U further increases up to 1/2 where we find a first- 
order transition from the ferromagnetic metal to the Mott 
insulator. 

Let us now turn to the case of vanishing Hund's rule 
coupling J = and U' = U. As shown in Fig. EH a quali- 
tatively different phase diagram is observed. For 77 = we 
saw in Sec. 14.11 that the paramagnetic metal is quite sta- 
ble due to the enhanced degeneracy of the lowest atomic 
configurations and that there is a joint transition to the 
Mott-insulating phase. This continuous transition is de- 
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noted by Pc in Fig. 1101 An arbitrarily small field 77 lowers 
the energy of the BI phase compared to the Mott insulator 
and therefore we find at the breakdown of the metallic so- 
lution for any finite 77 and U a first-order transition to the 
BI phase characterized by the mean field S2 = 1. Similar 
to a finite J, a finite 77 lifts the degeneracy of the six lowest 
on-site configurations, orbital fluctuations are suppressed 
and therefore the stability of the metallic phase is reduced 
with increasing crystal field. 

Note that our calculations simplify the true behavior of 
the system near the transition lines because the uniform 
mean-field approximation always reproduces the results 
of the atomic limit whenever the kinetic energy vanishes. 
Nevertheless they give some insight of the rich behavior of 
the system in the present of a crystal field which relatively 
shifts the atomic energy levels of the two orbitals. 



4.3.3 Mott transition away from half filling 

We now address the question of what happens away from 
half filling, 71 = 2 — 26, but again with zero crystal-field 
splitting. Particle-hole symmetry allows to concentrate on 
6 > 0. In general we observe a Mott transition in the nar- 
row band which lies at an increased interaction strength 
compared to the case 6 = (see Fig. lllfl . Because the sec- 
ond band is always away from half filling it stays metal- 
lic. As a representative example we show in Fig. El the 
band-renormalization factors Qa and the charges ria for 
fixed ratios U'/U = 1/2 and J/U = 1/4 and given dop- 
ing 6 = 0.03. Let us first look at the noninteracting case, 
U = 0. The ground-state energy per site in this case is 



Eg = eUl ~ Si) +e°2{l- Si) 



(17) 



where Sa is the deviation from half filling in band a and 
£° the average kinetic energy per site in band a for a 
half-filled band. Since e° is proportional to the band- 
width Da we find that the kinetic energy is optimized 
by choosing the charge imbalance A = (ni — n2)/2 = 
(Di - D2)/{Di + D2)S. For 6 = 0.03 and D1/D2 = 1/2 
this gives the value A — —0.01 as seen in Fig. El Thus, for 
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Fig. 12. Band-renormalization factor qa for finite hole doping 
5 = 0.03 and the charge distribution of the two bands for a 
fixed ratio U' /U = 1/2 and J/U = 1/4. 



U — 0, the narrow band serves as a charge reservoir that 
allows to bring the broader band closer to half filling. With 
increasing interaction the Coulomb energy causes a trans- 
fer of electrons from the wide band to the narrow band 
in order to reduce the intraorbital repulsion. Thus, with 
increasing interaction strength, the broader band serves 
as a charge reservoir. This gives rise to a half-filled band 
{A = (5) at a certain interaction strength and to an OSMT. 
The metallic behavior of the second band is due to the fi- 
nite hole doping and we find go = 2(5 up to first order in S 
in the atomic limit Sa/U ^0 p2 |. 



5 Discussion 

5.1 Comparison to known results at half filling 

On the mean-field level of the slave-boson approach we 
investigated the Mott transition in the two-band Hubbard 
model with different bandwidths and confirmed several 
aspects of previous work. As reported by several authors 
0,|^^^^^|li,|l|l3,|Ill we observe the OSMT 
and consequently an intermediate phase where only the 
narrow band is insulating whereas the wide band still has 
metallic properties. 

Our mean-field calculations predict that for U = U' + 
2 J the OSMI phase is limited to a small parameter regime 
in the U-U' phase diagram which is characterized by a 
rather high value of J (Fig. EJ. Compared to the phase 
diagram shown in i^] the strength of the OSMI phase is 
strongly reduced. In view of our treatment of the Hund's 
rule coupling this can be expected. As pointed out in 

pair- hopping and the spin-flip term of the full 
Hund's rule coupling lead to a stabilization of the OSMI 
phase. Since these terms are omitted in our calculations 
(and also in previous QMC studies 0,|l|lll) the OSMI 
phase is strongly reduced. Nevertheless, also with an Ising 
like Hund's rule coupling we can clearly resolve a sequence 
of individual Mott transitions in our slave-boson approach. 

In addition, we showed that on the mean-field level the 
PM-OSMI transition is second-order whereas the OSMI- 
MI and PM-MI transitions are first-order. Near these tran- 
sitions there coexist two different solutions and the en- 
ergy crossing results in a first-order transition (Fig. [Tj). 
The same behavior was reported in in the framework 



of the Gutzwiller approximation. Different results were 
found within other methods^ but non of the used meth- 
ods is rigorous and the order of the transitions remains an 
open problem. Furthermore, temperature as well as pair- 
hopping and spin-flip terms might affect the order of the 
phase transitions Q . 

For the case U = U' and J = orbital fluctuations lead 
to a stabilization of the metallic phase and for a fixed ratio 
of the bandwidths D1/D2 ~ 1/2 a joint second-order tran- 



sition is observed (Fig. 3(b) I. For the case of two bands 
of much different bandwidths, Di ^ D2, the Gutzwiller 
approximation and related mean-field theories predict the 
existence of an OSMI if the ratio D1/D2 is below a critical 
value. This was first reported in DMFT calculations 

give clear evidence that the localized band is not a con- 
ventional Mott insulator but has spectral weight down to 
arbitrarily small energies. However, an instability toward 
an orbitally ordered phase might play an important role 
for U = U' and should also be taken into account in future 
investigations. 

5.2 Results for shifted bands 

In real materials different atomic orbitals are usually not 
degenerate so that each band has non-commensurate fill- 
ing. This extension of our model has lead to a considerably 
richer phase diagram. Such models correspond to the sit- 
uation found in Ca2-a;Sra;Ru04 which has three bands of 
partial filling. This material has been an initial motivation 
for the study of the OSMT Il3|. Anisimov and coworkers 
proposed that the Ca-Sr substitution varies band param- 
eters which in the end leads to a Mott transition in two 
of the three bands ^ further example which belongs 
likely to this class is the compound FeSi. This compound 
is a small gap semiconductor 0|. On the other hand, 
replacing Ge for Si gives rise to a ferromagnetic metal. 
Alloying FeSii-^^Ge^^ allows in principle for a continuous 
change of the band parameters such that the transition 
can be observed. However, the transition is simultaneous 
accompanied by an abrupt transition in the crystal lattice 

Within the mean field approximation we find the fol- 
lowing situation. For small crystal field splitting an OSMT 
is observed. In general, the Mott transition is shifted to 
higher values of the interaction parameters. Due to the 
crystal field, a band insulating phase and, in the presence 
of the Hund's rule coupling, also a ferromagnetic phase ap- 
pear in the phase diagram (Fig. ^ . In the ferromagnetic 
phase, a few electrons populate the upper band with a fi- 
nite net magnetization. For the case U = U' and J = we 
find a totally different behavior fFig.ll()|l. With increasing 
field, the metallic phase is less stable because the crystal 
field suppresses orbital fluctuations, similar to a flnite J, 
by breaking the degeneracy of the local states. 

For flnite doping our calculations suggest that there 
is in general a Mott transition in the narrow band for 
not too strong doping. Although strongly correlated, the 
second band stays metallic due to the flnite doping. This 
was also reported in 
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5.3 Conclusions 

In summary, the mean-field theory based on the slave- 
boson approach of Kotliar and Ruckenstein gives results 
which are in good qualitative agreement with DMFT cal- 
culations. While we restricted ourselves to density-density 
interactions, our discussion provides a simple physical pic- 
ture in a wide range of parameters. The transverse spin 
coupling and the on-site interorbital pair hopping had to 
be dropped for practical reasons. Nevertheless, we believe 
that the effects are rather of quantitative than qualitative 
nature. 

The method used emphasizes the on-site correlation 
and intersite correlations remain treated at a minor level 
only. Thus we have ignored symmetry breaking instabil- 
ities which double the size of the unit cell, such as anti- 
ferromagnetic instabilities or orbital order. These orders 
depend strongly on the detailed band structures and cou- 
pling topologies. In most of our discussion, however, we 
neglected the band structure aspect. Obviously nesting 
properties would play a major role in this context. Generic 
bands without nesting, however, follow more likely the 
"plain" behavior of the simple flat density of states mod- 
els that we discussed here. It would be interesting in fu- 
ture studies to extend the scheme by including also band 
structure effects and the related ordering phenomena. 
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Leuenberger, A. Georges and M. Ferrero for stimulating discus- 
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